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ABSTRACT 

We summarize some of the difficulties that confront lattice calculations of 
non-leptonic kaon decay matrix elements. We review some of the methods that 
have been proposed to overcome these difficulties, and discuss the importance 
of one-loop ChPT in this respect, including the role of 0(p 4 ) operators. 



In this talk, we will discuss some aspects of lattice computations of 
hadronic matrix elements, in particular those relevant for non-leptonic kaon 
decays. A more extensive review, with additional references, can be found in 
ref. jj]. In general, Lattice QCD provides us with euclidean correlation func- 
tions which are typically computed at unphysically large values of quark masses 



speaker at conference 



and/or unphysical momenta (see below). In order to extract physical informa- 
tion, we can analytically continue the unphysical lattice results using chiral 
perturbation theory (ChPT), at least when we are concerned with the physics 
of Goldstone bosons of QCD. 

There are several methods which are currently being used or investigated 
for the computation of kaon-decay rates. The oldest is to compute the K tttt 
matrix elements with not only the kaon but also the pions at rest || , and use 
ChPT in order to continue to physical (i.e. energy-conserving) momenta. Al- 
ternatively, K —> 7T7T amplitudes can also be related through soft-pion theorems 
to the simpler K — > 7r matrix elements which are, at least in principle, easier to 
compute on the lattice || (see also ref. Q) Recent numerical work employing 
both these methods can be found in refs. || ||. More recently, it has been 
suggested that the finite size of the spatial volume can be used as a tool to get 
to the physical K — > tttt matrix elements [fr| . 

In what follows we will always assume that we are dealing with continuum- 
extrapolated lattice results, so that the usual ChPT techniques directly apply. 
On the lattice, one extracts K — > tttt matrix elements from euclidean correlation 
functions of the form 

C(t 2 ,h) = (0|7r(g,t 2 )7r(-g,t 2 )O weafe (t 1 )i<r(0)|0) (1) 
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for t\ large. The dominant contribution to this correlation function at large 
euclidean times comes from the state with both pions at rest, i.e. q = 
and Eg = 2m w + 0(1/L 3 ), while the desired physical state, which has \q\ = 
| y/ m 2 K — 4m^ , is buried in the tower of excited states of eq. (0) || . Note that 
for \n) — |7r(<f = 0)7r(g = 0)) energy is injected by the weak operator. 

These simple observations lead us to the three aforementioned meth- 
ods. In the approach proposed in ref. one may work in a finite volume, 
in which the smallest non-zero momentum \q\ — 2it / L is chosen such that 
rriK = Ei ~ 2 \J m% + q 2 (omitting interaction effects). If one can single out 
the first excited state, this gives the desired matrix element at finite volume, to 
which a correction factor can be applied to convert to infinite volume 0. At 
the physical values of meson masses, this requires a rather large volume of order 
m^L w 4, or L « 6 fm. It will be interesting to see whether this method can 
be made to work in practice. This proposal together with a related approach 



M is further discussed by Guido Martinelli in these proceedings || . 

A second, technically much simpler, method is to compute the unphysical 
matrix element (Tr(q — 0)7r(<f — Q)\O wea k\K), and use ChPT to convert it to 
the physical one Q]. This approach overcomes the problems implied by the 
Maiani- Testa theorem ]8j since in the unphysical configuration with both pions 
at rest no final-state interaction phase is generated. A third method is based 
on the observation that ChPT relates K — + irir matrix elements to K — > 7r and 
K — > transitions || . Issues which arise in the use of the latter two methods 
are: 1) the size of chiral corrections, and related, 2) the role of 0(p A ) ChPT 
operators, 3) finite-volume effects, and 4) quenching. Issues 3 and 4 can also 
be investigated through the use of one-loop ChPT. 

As an example, let us consider the A J = 3/2 decay K + — ► 7r + 7r°. To one 
loop the physical matrix element is given by 



[K+ -> Tt + ir\ hys 



rn 



K 



1 + chiral logs + dx + dn 
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while for the unphysical matrix element with both pions at rest and Mr 
M n — M the degenerate meson mass on the lattice jlO) one has 

N M 2 + M, 



(2) 



[K + -> ir + ir°) unph y S ^ a 27 (7V) 



M 2 



1 + chiral logs + d s 



sea 
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N is the number of (degenerate) dynamical fcrmions (=sea quarks) on the lat- 
tice, with corresponding Goldstone-boson mass M sea (which on the lattice does 
not have to be equal to the valence- meson mass M). The quenched approx- 
imation corresponds to N = (see rcf. jl| for more details). The low-energy 
constant (LEC) a 2 '' (N) is that of the real world for TV = 3, even if the quark 
masses on the lattice are not at their physical values [0 . In the fully quenched 
theory, it may have a different value. Similar considerations also apply to the 
linear combinations of 0(p A ) LECs dK,n,s,v In addition, the linear combina- 
tions showing up in the unphysical matrix element can be different from those 
of the physical one: For arbitrary (on-shell) momenta px , Pm and Ptt 2 the 
contribution from 0(p 4 ) operators looks like 

(tree level) x — — ^ (AM 2 K + BM 2 + Cp n ■ p„ 2 + Dp K ■ (p ni + p X2 )) . (4) 



In the physical case px — V-kx + V-ki so that the expression in parentheses 
becomes (A — ~C — D)M K + (B + C)M 2 , whereas in the unphysical case 
p lxi = pir 2 = (0, iM n ) and px — (0, iMk) lead to the combination AM K + (B — 
C)M 2 - 2DM K M^, which does not yield the same information. The upshot 
is that one can determine the leading-order LECs, such as a 27 , but not easily 
all needed 0(p 4 ) LECs.|] The "chiral logarithms" in eqs. (0j3) are different in 
the two cases, and, for typical lattice quark masses, they are sizable. This is 
an indication that one should at least use ChPT to one loop in the analysis of 
lattice results for K — > tttt matrix elements. 

To conclude this example, we note that there are also power-like finite- 
volume corrections, here given for a spatial box L 3 with periodic boundary 
conditions. These corrections come from pion-rescattering diagrams They 
may be large: for / = 160 MeV, M n = 500 MeV and M W L = 6 (typical lattice 
values of these parameters) , they are about 20% of the tree- level contribution. 
By comparing computations on different volumes, JLQCD found good agree- 
ment between the prediction from ChPT of cq. (||) and results from lattice 
QCD (first paper of Q). In general, one- loop ChPT gives a good fit of the 
JLQCD results (which are quenched, at a -1 = 2 GeV), with a reasonable value 
of d v . It appears that one-loop ChPT is sufficient to explain the discrepancy 
[jl2| between unphysical (K — > Trir)Ai=3/2 lattice results converted to physical 
ones at tree level only of ChPT and the experiment. 

A similar analysis can be carried out for AI = i decays. This case is 
however more complicated for a variety of reasons (0(p A ) LECs being only one 
of them - for a nice review see ref. JH|), and it may be advantageous to use the 
chiral-symmetry connection to K — > 7r matrix elements. At tree level in ChPT, 
this works as follows 0. The desired K — > nn matrix elements are given by 

4,- 

(TT + n-\27plet\K°) = ( m l- m l)a 27 , (5) 

r 
Ai 

{-K + Ti-\octet\K a ) = - — { m 2 K -ml)a\ , 

in terms of the 0(p 2 ) (27^,1^) and (8l,1r) LECs a 27 and af respectively. 

1 In the more complicated case of AI = | decays, also total-derivative 
operators, such as <9 A1 [£9 p £' 1 ', EM* ± MT,^]±,is, appear at 0(p 4 ). 



They can also be obtained from 
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(Tr + \27plet\K+) = -—M 2 a 27 , (6) 
J 

{■n + \octet\K + ) = j^M 2 (al - a%) , 

taking Mk — M^ = M on the lattice. The new (unphysical) LEC a\ can be 
determined from the matrix element 

(0\octet\K°) = -{M\ - M 2 )a 8 2 . (7) 

To lowest order in ChPT, K — » tttt matrix elements can be therefore obtained 
from K — > ir and K — > 0, but the likely size of chiral corrections on the lattice 
(and in nature) makes it imperative to extend the analysis to, at least, next- 
to-leading order in ChPT. At one loop, ChPT predicts a behavior like JT3] 

4M 2 

(ir+\octet\K+) = — [a?(l + logs) - ag(l + logs) + A V M 2 + A s NM 2 sea ] , 

(8) 

and analogous expressions for the other matrix elements. Indeed, for typical 
lattice masses, the logarithms can be large. Hence, they and the polynomial 
terms with coefficients A v . s , which are linear combinations of 0(p A ) LECs, 
will have to be taken into account in fits to lattice data. Unfortunately, the 
combinations of 0(p 4 ) LECs contained in K — > tt matrix elements do not carry 
enough information to fully determine physical K — > 7T7t matrix elements up to 
0{p A ). Also note that for K — > and K — > tt with Mk ^ M^ total-derivative 
operators can (and do) contribute. In this context we mention that there also 
exist subtleties with penguin operators in (partially) quenched QCD which 
have been discussed in ref. |jT5f . 

One lesson to be learned is that, in our view, it should be possible to 
extract reliable values of the 0(p 2 ) LECs a\ and a 27 from the lattice computa- 
tion of K — > 7r transition matrix elements. This is interesting in its own right, 
specially since phenomenological estimates of these LECs already exist . In 
addition, we believe that the analysis of lattice results within the approaches 
considered here requires the use of next-to-leading order (i.e. 0(p A )) ChPT for 
a variety of reasons. The most important of these are to check for convergence 
of the chiral expansion, to reduce systematic errors coming from the use of 
ChPT, and to understand finite-volume effects. 
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